DIFFEOMORPHIC APPROXIMATION OF 
SOBOLEV HOMEOMORPHISMS 



TADEUSZ IWANIEC, LEONID V. KOVALEV, AND JANI ONNINEN 

Abstract. Every homeomorphism h: X — > Y between planar open sets 
that belongs to the Sobolev class W 1,P (K, Y), 1 < p < oo, can be ap- 
proximated in the Sobolev norm by "if °°-smooth diffeomorphisms. 



1. Introduction 

By the very definition, the Sobolev space W 1,P (1L, M), 1 ^ p < oo, in 
a domain X C M n , is the completion of -smooth real functions having 
finite Sobolev norm 

INI^^x) = ll n llifp(X) + ||Vu||jsf P (x) < oo. 
The question of smooth approximation becomes more intricate for Sobolev 
mappings, whose target is not a linear space, say a smooth manifold [HJdHl 
[201 [21] or even for mappings between open subsets X, Y of the Euclidean 
space M n . If a given homeomorphism h: X Y is in the Sobolev class 
W 1,P (K, Y) it is not obvious at all as to whether one can preserve injectivity 
property of the ^^-smooth approximating mappings. It is rather surpris- 
ing that this question remained unanswered after the global invertibility of 
Sobolev mappings became an issue in nonlinear elasticity [U |T71 [3TJ [35] . It 
was formulated and promoted by John M. Ball in the following form. 

Question. [6l[7] If h £ 'W 1,P (K, W 1 ) is invertible, can h be approximated in 
W ,p by piecewise affine invertible mappings? 

J. Ball attributes this question to L.C. Evans and points out its relevance 
to the regularity of minimizers of neohookean energy functionals [5j El EJ 
[16l 134] . Partial results toward the Ball- Evans problem were obtained in |30j 
(for planar bi-Sobolev mappings that are smooth outside of a finite set) and 
in [TO] (for planar bi-Hdlder mappings, with approximation in the Holder 
norm) . The articles (6J [33] illustrate the difficulty of preserving invertibility 
in the approximation process. In [21] we provided an affirmative answer 
to the Ball- Evans question in the planar case when p = 2. In the present 
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paper we extend the result of [M] to all Sobolev classes W 1,P (K,Y) with 
1 < p < oo. The case p = 1 still remains open. 

Let X be a nonempty open set in IR 2 . We study complex- valued functions 
h = u + iv : X -> C ~ R 2 of Sobolev class W^ P {X, C), 1 < p < oo. Their 
real and imaginary part have well defined gradient in J£ P (X,R 2 ) 

Vu : X -»• R 2 and : X -»• M 2 . 

Then we introduce the gradient mapping of /i, by setting 

(1.1) Vh = (Vu,Vv): X R 2 x M 2 . 

The jSf p -norm of the gradient mapping and the p-energy of h are defined by 



(1.2) \\Vh\\^ P{x) 



(\Vu\ p + \Vv\ p ) 



E x [h] = E p [h] = ||Vfc" p 



The reader may wish to notice that this norm is slightly different from what 
can be found in other texts in which the authors use the differential matrix 
of h instead of the gradient mapping, so 



(i.3) \\Dh\y P{x) 



(|Vit| 2 + \Vv 



p 



Thus our approach involves coordinate-wise p-harmonic mappings, which we 
still call p-harmonic for the sake of brevity. We shall take an advantage of 
the gradient mapping on numerous occasions, by exploring the associated 
uncoupled system of real p-harmonic equations for mappings with smallest 
p-energy. Our theorem reads as follows. 

Theorem 1.1. Let h: X 22*2. y be an orientation-preserving homeomor- 
phism in the Sobolev space W^f(X,Y), 1 < p < oo, defined for open sets 
X,¥cR 2 . Then there exist tf 00 -diffeomorphisms h e : X ^ Y, £ = 1, 2, . . . 
such that 

(i) h e -h£ W hp (X, R 2 ), £ = 1,2,... 

(ii) lim (hg — h) = 0, uniformly on X 

>-oo 

(Hi) lim \\Vhn — Vh\\^p(X) = 

(iv) \\Vh£\\jz> P( x) «S ||V/i||^ P(x) , for £ = 1,2,... 

(v) If h is a c £ co -diffeomorphism outside of a compact subset of X, then 
there is a compact subset o/X outside which hi = h, for all £ = 1,2,... 

A straightforward triangulation argument yields the following corollary. 

Corollary 1.2. Let h: X 2H!$ Y be an orientation-preserving homeomor- 
phism in the Sobolev space J0£*(X,Y), 1 < p < oo, defined for open sets 
X,¥cR 2 . Then there exist piecewise affine homeomorphisms hg: X Y, 
£ = 1,2,... such that 

(i) h e -heWo' p (X,R 2 ), £ = 1,2,... 
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(ii) lim (he — h) = 0, uniformly on X 

i— >oo 

(Hi) lim II Vht — V/ill^fpfxi = 0. 

(raj ///i is affine outside of a compact subset o/X ; i/jen there is a compact 
subset o/X outside which he = h, for all I = 1,2,... 

We conclude this introduction with a sketch of the proof. The construc- 
tion of an approximating diffeomorphism involves five consecutive modifi- 
cations of h. Steps 1, 2, and 4 are p-harmonic replacements based on the 
Alessandrini-Sigalotti extension [3] of the Rado-Kneser-Choquet Theorem. 
The other steps involve an explicit smoothing procedure along crosscuts. 
For this, we adopted some lines of arguments used in J. Munkres' work [32j . 

2. p-HARMONIC MAPPINGS AND PRELIMINARIES 

Let O be a bounded domain in the complex plain C ~ M 2 . A function 
u: Q, — > K in the Sobolev class 1 < p < oo, is called p-harmonic if 

(2.1) div \Vu\ p ~ 2 Vu = 
meaning that 

(2.2) f (\Vu\ p - 2 Vu,Vip) = for every ip £ ^(fi). 

The first observation is that the gradient map / = Vm: $7 — > M 2 is K- 
quasiregular with 1 ^ K ^ max-fj? — 1,1/ (jp — 1)}, see |12j . Consequently 
u G ^ip) with some < a = a(p) ^ 1. In fact [25] the foremost 
regularity of a p-harmonic function (p ^ 2) is ^^"(O), where the integer 
k ^ 1 and the Holder exponent a S (0, 1] are determined by the equation 

7p - 6 + Jp 2 + Yip - 12 1 1 

k + a = > 1 -I — . 

6p — 6 3 

Thus, regardless of the exponent p, we have u S ^ lo '"(^) with a = 1/3. 
Clearly, by elliptic regularity theory, outside the singular set 

S = {z€n-. Vu(z) = 0}, 

we have u £ < ^'° o (0,\S). The singular set, being the set of zeros of a quasireg- 
ular mapping, consists of isolated points; unless u = const. Pertaining to 
regularity up to the boundary, we consider a domain Q whose boundary 
near a point z £ Oil. is a ^^-smooth arc, say T C d£l. Precisely, we as- 
sume that there exist a disk D = D(z Q , e) and a ^^-smooth diffeomorphism 
ip : D C such that 

<p(D n O) = C+ = {z : Im z > 0} 

«p(r) = M = {z: Imz = 0} 
ip(D \ O) = C_ = {z : Im z < 0}. 
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Proposition 2.1 (Boundary Regularity). Suppose u G W 1,p {^l) n ^(Q) is 
p-harmonic in Q and c to°° -smooth when restricted to T. Then u is c £ l ' a - 
regular up to V , meaning that u extends to D as a 'tf > a (D) -regular function, 
where a depends only on p. 

2.1. The Dirichlet problem. There are two formulations of the Dirichlet 
boundary value problem for p-harmonic equation; both are essential for our 
investigation. We begin with the variational formulation. 

Lemma 2.2. Let u Q G W l > p {Vt) be a given Dirichlet data. There exists 
precisely one function u G u Q + W which minimizes the p-harmonic 

energy: 

£ p [u] = inf I J \Vw\ p : w G u Q + #^(ft)J . 

The solution u is certainly a p-harmonic function, so ^^'"(il) -regular. 
However, more efficient to us will be the following classical formulation of 
the Dirichlet problem. 

Problem 2.3. Given u G ^(dCl) find a p-harmonic function u in $7 which 
extends continuously to O such that m gQ = u a . 

It is not difficult to see that such solution (if exists) is unique. How- 
ever, the existence poses rather delicate conditions on d£l and the data 
u G ^(fl). We shall confine ourselves to Jordan domains J! C C and the 
Dirichlet data u Q G ^(O) of finite p-harmonic energy. In this case both for- 
mulations are valid and lead to the same solution. Indeed, the variational 
solution is continuous up to the boundary because each boundary point of a 
planar Jordan domain is a regular point for the p-Laplace operator A p \18\ 
p. 418]. See [22l 6.16] for the discussion of boundary regularity and relevant 
capacities and \27\ Lemma 2] for a capacity estimate that applies to simply 
connected domains. 

Proposition 2.4 (Existence). Let f2 C C be a bounded Jordan domain 
and u G W 1,p (Sl) n <Jf (H). There exists, unique, p-harmonic function u G 
W x v{Vl) n ^(0) such that u ]an = u o]gn . 

2.2. Rado-Kneser-Choquet Theorem. Let h = u + iv be a complex har- 
monic mapping in a Jordan domain U that is continuous on U. Assume that 
the boundary mapping h : <9U T is an orientation-preserving homeomor- 
phism onto a convex Jordan curve. Then h is a ^°°-smooth diffeomorphism 
of U onto the bounded component of C \ V. Thus, in particular, the Ja- 
cobian determinant J(z,h) = \h z \ 2 — |/i s | 2 is strictly positive in U, see [T3J 
p. 20]. Suppose, in addition, that <9U contains a ^^-smooth arc 7 C 9U, 
and h takes 7 onto a ^"-smooth subarc in T. Then h is ^^-smooth up 
to 7 and its Jacobian determinant is positive on 7 as well, see [151 p. 116]. 
Numerous presentations of the proof of Rado-Kneser-Choquet Theorem can 
be found, [15]. The idea that goes back to Kneser [26] and Choquet [13] 
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is to look at the structure of the level curves of the coordinate functions 
u = Re/i, v = Imh and their linear combinations. These ideas have been 
applied to more general linear and nonlinear elliptic systems of PDEs in 
the complex plane [H] , see also [U 0, ESI EH] for related problems concerning 
critical points. In the present paper we shall explore a result due to G. 
Alessandrini and M. Sigalotti j3] for a nonlinear system that consists of two 
p-harmonic equations 



Call it uncoupled p-harmonic system. The novelty and key element in [3] is 
the associated single linear elliptic PDE of divergence type (with variable co- 
efficients) for a linear combination of u and v. Such combination represents a 
real part of a quasiregular mapping and, therefore, admits only isolated criti- 
cal points. We shall not go into their arguments in detail, but instead extract 
the following p-harmonic analogue of the Rado-Kneser-Choquet Theorem. 

Theorem 2.5 (G. Alessandrini and M. Sigalotti). Let U be a bounded Jor- 
dan domain and h = u + iv.U — > C be a continuous mapping whose coor- 
dinate functions u, v € W 1,P (U), 1 < p < oo, are p-harmonic. Suppose that 
h: <9U ^2-4 7 is an orientation-preserving homeomorphism onto a convex 
Jordan curve 7. Then 

(i) h is a c rf°° -diffeomorphism from U onto the bounded component o/C\7. 
In particular, 



(ii) If, in addition, <9U contains a 'it? 00 -smooth arc V C <9U and h(T) is 
a ff 00 -smooth subarc in 7, then h is ( ta 1,a -regular up to T, for some 
< a = a(p) < 1 (actually c € ca ). Moreover J(z, h) > onT as well. 

This theorem is a straightforward corollary of Theorem 5.1 in [3J. How- 
ever, three remarks are in order. 

(1) In their Theorem 5.1 the authors of [3] assume that U satisfies an 
exterior cone condition. This is needed only insofar as to ensure 
the existence of a continuous extension of a given homeomorphism 
$ : d\J — > 7 into U whose coordinate functions are p-harmonic in 
U. Obviously, such an extension is unique, though the p-harmonic 
energy need not be finite. Once we have such a mapping the exterior 
cone condition on U for the conclusion of Theorem 5.1 is redundant, 
see Remark 3.2 in [3]. This is exactly the case we are dealing with 
in Theorem 12.51 

(2) In regard to the statement (ii) we point out that in Theorem 5.1 
of [3J the authors work with the mappings that are smooth up to the 
entire boundary of U. Nonetheless their proof that J(z, h) > on 
dU is local, so applies without any change to our case (ii). 




1 < p < 00, h = u + iv. 



J(z,h) = \h z \ 2 - \h s \ 2 > (till 
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(3) Since J(z, h) > in U up to the arc T C d\J the coordinate functions 
of h have nonvanishing gradient. This means that p- harmonic equa- 
tion is uniformly elliptic up to T. Consequently, h is ^°°-smooth on 
U up to r. 

2.3. The p-harmonic replacement. Let Q, be a bounded domain in M 2 ~ 
C. We consider a class A(fl) = A p (£l), 1 < p < oo, of uniformly continuous 
functions h = u + iv: Q — >• C having finite p-harmonic energy and furnish it 
with the norm 

INUp(q) = IWk(fi) + \\^ h \\j?p(n)- 
The closure of ^ o °°(0) in A p {n) will be denoted by A%(n). 

Proposition 2.6. Let U (s SI be a Jordan subdomain of £1. There exists a 
unique operator 

R w : A p (n) -^A p (tt) 
(nonlinear if p ^2) such that for every h £ A p (£l) 

Ry/i = h in O \ U 

(2.3) Ru G h + Wj-' P (V) 

ApRu/i = inU 

(2.4) E n [Rufc] < E n [h] 



Equality occurs in (2.4) if and only if h is p-harmonic in U. 
Proof. For h = u + iv we define 

Rnj/i = Ruw + i Ru^- 

It is therefore enough to construct the replacement for real- valued functions. 
For u G ^l p (0) real, we define 



u in \ U 
u in U 

where u is determined uniquely as a solution to the Dirichlet problem 

fdiv|Vu| p ~ 2 Vn = in U 



so conditions (2.3) are fulfilled. That Run is continuous in Q is guaranteed 



by Proposition 2.4 The solution u is found as the minimizer of the p- 
harmonic energy in the class u + Wo' p (\3), so we certainly have 

En [Run] ^ E n [u} 

The same estimate holds for the imaginary part of h, so adding them up 
yields 

En \l;h < E n [h]. □ 

Remark 2.7. The reader may wish to know that the operator Ru: A(£l) — 
A(Cl) is continuous, though we do not appeal to this fact. 
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2.4. Smoothing along a crosscut. Consider a bounded Jordan domain 
U and a ^^-smooth crosscut T C U with two distinct end-points in <9U. By 
definition, this means that there is a ^^-diffeomorphism ip: C U such 
that r = </?(M), and its distinct endpoints are given by 

lim <p(x) G dU 

x—t—ao 

lim (fix) 6 9U 

x— >oo 

Such r splits U into two Jordan subdomains 

U+ = ^(C+), C + = {z: ImoO} 
U_ = 99(C_), C_ = {z: Imz < 0}. 

Suppose we are given a homeomorphism /: U — > C such that each of two 
mappings 

/ : U+ ->• R 2 and / : U_ ->• M 2 
is ^"-smooth up to T. Assume that for some constant < m < oo we have 

\Df(z)\ < m and detD/(z) ^ — 

m 

on U+ and on U_. Thus / : U — > M 2 is in fact locally bi-Lipschitz. 

Proposition 2.8. Under the above conditions there is a constant < M < 
oo such that for every open set V C U containing T one can find a home- 
omorphism g: U /(U) which is a < ^°° -diffeomorphism in U, with the 
following properties: 

(2.5) g[z) = f(z), forze (U\V)UT 

(2.6) \Dg(z)\ ^ M and detDg(z) > — on U. 

The key element of this smoothing device is that the constant M is inde- 
pendent of the neighborhood V of T, see Figure [TJ The proof is given in |24j 
following the ideas of |32j. 

We shall recall similar smoothing device for cuts along Jordan curves. 
Let U be a simply connected domain with ^^-regular cut along a Jordan 
curve r C U. This means there is a diffeomorphism ip: C U such that 
r = ifiS 1 ), S 1 = {z £ C: \z\ = 1}. As before T splits U into 

V+=ip(B + ), D + = {z: \z\ < 1} 
U_=^(D_), B_ = {z: \z\ > 1}. 

Suppose we are given a homeomorphism / : U — > M? such that each of two 
mappings 

/:U+^R 2 and /:U_^R 2 
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Figure 1. Jordan domain with a crosscut T and its neigh- 
borhood V. 



is 'rf^-smooth up to T. Assume that for some constant < m < oo we have 

1 



|.D/(2)| < m and det Df(z) ^ 



on U+ and U_ . 



Proposition 2.9. Under the above conditions there is a constant < M < 
oo such that for every open set V C U containing T one can find a c to°°- 
diffeomorphism g: U /(U) with the following properties 

(2.7) g(z) = f(z),forze(U\Y)ur 

(2.8) \Dg(z)\ < M and det Dg(z) > ^ on U. 

Having disposed of the above preliminaries we shall now proceed to the 
construction of the approximating sequence of diffeomorphisms. 

3. The proof 

3.1. Scheme of the proof. Let us begin with a convention. We will often 
suppress the explicit dependence on the Sobolev exponent 1 < p < oo in 
the notation, whenever it becomes selfexplanatory. For every e > we shall 
construct a ^°°-diffeomorphism h: X Y such that 

(A) H-h€Ao(X) 

(B) \\h-h\\ VQ Q 

(c) \\vh-vh\y P(x) ^ e 

(D) E x [ft] < E x [h] 
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(E) If h is a ^°°-diffeomorphism outside of a compact subset of X, then 
there exist a compact subset of X outside of which we have h = h, for 
all e > 0. 

We may and do assume that h is not a ^°°-diffeomorphism, since other- 
wise h = h satisfies the desired properties. Let x Q £ X be a point such that 
h fails to be ^°°-diffeomorphism in any neighborhood of x . 

We shall consider dyadic squares in Y with respect to a selected rectangu- 
lar coordinate system in M? . By choosing the origin of the system we ensure 
that h(x Q ) does not lie on the boundary of any dyadic square. 

Let us fix e > 0. The construction of h proceeds in 5 steps, each of 
which gives a homeomorphism h~ k : X Y, k = 0, 1, . . . , 5, in the Sobolev 
class y^f(X, Y) such that Hq = h, % G h k -i + „4 (X), k = 1,...,5 and 
^5 = h is the desired diffeomorphism. For each k = 1, 2, . . . , 5 we will secure 
conditions analogous to (A)-(E). Namely, 
(A k ) H k - %_i G Ao(X) 
(B k ) \\H k - %_i||^(x) ^ e/5 
(C fc ) ||V^fc - Vft fc -i||jS!*(X) < e/5 
(D k ) || V^i||j£fp(x) < ||V^o||^fp(x) - 25, for some 5 > 0; 

||V%||^f P (x) «S ||V%_i||^p( X ), for k = 2,4; 

l|V%||^ P( x) «S ||V%_i||^.p(x) +5, for fc = 3,5 
(-Efe) If h k -\ is a ^^-diffeomorphism outside of a compact subset of X, then 

there exists a compact subset in X outside which we have fi k = h k -i 

for all e > 0. 

3.2. Partition of X into cells. Let us distinguish one particular Whitney 
type partition of Y and keep it fixed for the rest of our arguments. 

oo 

Y=\J%, 

v=l 

where Y^ are mutually disjoint open dyadic squares such that 

diamY,, < dist(Y„, dY) sC 3diamY^ for v = 1, 2, . . . 

unless Y = R 2 , in which case Y^ are unit squares. Thus the cover of Y by 
Y u is locally finite. The preimages 

X v = h-\Y U ), e/ = 1,2, ... 

are Jordan domains which we call cells in X. In the forthcoming Step 1 we 
shall need to further divide each cell into a finite number of daughter cells 
in X. Note that all but finite number of cells X u , v = 1,2, ... lie outside a 
given compact subset of X. 

Step 1 

To avoid undue indexing in the forthcoming division of cells, we shall 
argue in two substeps. 
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Step la. Examine one of the cells in X, say X = X^, for some fixed v = 
1,2,.... Call it a parent cell. Thus h(X) = T is the corresponding Whitney 
square T = ¥ v C ¥. To every n = 1, 2, . . . , there corresponds a partition of 
T into 4 n -dyadic congruent squares T j , i = 1 , . . . , 4 n 



T = Ti U • • • U T 4 n. 
This gives rise to a division of X into daughter cells Xj = h~ l (Ti 

x = XT u x~2 u • • • u x~4^. 

We look at the homeomorphisms 

/i:X~^T~, i = l,2,...4 n 



By virtue of Proposition 2.6 we may replace them with p-harmonic homeo- 
morphisms 

h i = R x .h:X i ^T~ i , i = l,2,...,A n 
which coincide with h on dXi. This procedure may not be necessary if 
h: 3ii — > Ti is already a ^"^-diffeomorphism. In such cases we always use 
the trivial replacement hi = h. After all such replacements are made, we 
arrive at a homeomorphism 

which is a ^^-diffeomorphism in each cell Xj and coincides with h on dXi. 
Obviously, 

h = h + Y^fti -h] £h + A a (X) 

2=1 

where [hi — h] Q stands for zero extension of hi — h outside X 4 - and, therefore, 
belongs to A Q (Xi). Furthermore, by principle of minimal p-harmonic energy, 
we have 

41 4" 

E x [h] = £ e Xi n < Yl = W 

i=l i=l 

The eventual aim is to fix the number of daughter cells in X. For this we 
vary n and look closely at the resulting homeomorphisms, denoted by f n . 
This sequence of mappings is bounded in A(X). It actually converges to 
h uniformly on X. Indeed, given any point x G X, say x G Xj, for some 
z = 1,2,..., 4 n , we have 

|/n(^) — h(x)\ = \hi(x) — h(x)\ ^ diamTi = 2~ n diamY. 

Thus 

lim f n = h, uniformly in X. 

On the other hand the mappings /„ are bounded in the Sobolev space 
W 1,P (X), so converge to h weakly in W l,p {X). The key observation now 
is that 

l|V/i||j^ ( x) < liminf HV/nll^aq < \\^h\\^ P{x) 
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because of convexity of the energy functional. This gives 
J^llV/nll^x) = \\Vh\\j? P{x) 

Then, the usual application of Clarkson's inequalities in Jzf p -spaces, 1 < p < 
oo, yields 

lim ||V/„ - Vh\y p(x) = 

meaning that f n — h — > in the norm topology of A(3C). We can now 
determine the number n = n u = n(X), simply requiring the division of X be 
fine enough to satisfy two conditions. 

fdiamTi = 2-™diamT ^ e/5, i = l,...,4 n 



(3.1) 



l|V/„ - Vh\\xp m < 5^7 



where we recall that 3t stands for X^, . 

Step lb. Now, having n = n u fixed for each cell 3C U , we construct our first 
approximating mapping 

h x : X 22$ Y 

by setting 

oo 

h 1 -.= h + J2lfn„ -h] Q eh + A (X) 

u=l 

where, as always, [/ n „ — h] stands for the zero extension of f Uv — h outside 
H u . This mapping is a ^°°-diffeomorphism in every daughter cell. Clearly, 
we have the condition 

(Ai) Hi — he A>(X). 



5 



oo oo 1 



Moreover, by the condition in (3.1) imposed on every n v , 
(Bi) - /»||»(x) < sup {diamTV Ti C Y„,i = 1,. ..,4 n "} < J 

and 

(C X ) IIV* - Vfc||^ w = - ™UW,) < (i/ U 2^ < U 

I/=l I/=l 

Regarding condition (-Di), we observe that summing up the energies over 
all daughter cells j£j C X^, i = 1, 2, . . . 4 n " and u = 1,2, . . . , gives the total 
energy of not larger than that of /i. Even more, since h fails to be a < ^°°- 
diffeomorphism in at least one of these cells, the p-harmonic replacement 
takes place in this cell and, consequently, h\ has strictly smaller energy. 
Hence 

(Z?i) ||V%||j^p(x) ^ ||V/t||j^ P (x) — 25, for some 5 > 0. 

Regarding condition (Ei), we note that under the assumption therein we 
made only a finite number of nontrivial (p-harmonic) replacements. The 
same remark will apply to the subsequent steps and will not be mentioned 
again. The step 1 is complete. 



12 TADEUSZ IWANIEC, LEONID V. KOVALEV, AND JANI ONNINEN 

Before proceeding to Step 2, let us put all daughter cells in X in a single 
sequence 

X , X , * * * CZ X. 

Thus from now on the daughter cells from different parents are indistin- 
guishable as far as the mapping fi\ is concerned. The point is that fi\ is a 
^^-diffeomorphism in every such cell, a property that will be pertinent to 
all new cells coming later either by splitting or merging the existing cells. 
Note that the images T Q = h(X a ), a = 1,2, ... , form a partition of Y into 
dyadic squares 

oo 

Y = I IT", where diamT" sC -. 
, 5 

a=l 




Figure 2. hi is a ^^-diffeomorphism in each cell X a C X. 

Step 2 

Step 2a. (Adjacent cells) Let C(Y) C Y be the collection of all corners 
of dyadic squares T a , a = 1,2,..., and V(X) C X denote the set of their 
preimages under h, called vertices of cells. Whenever two closed cells X a 
and X@, a / (3, intersect, their common part is either a point in V(X) or 
an edge, that is, a closed Jordan arc with endpoints in V(X). In this latter 
case we say that X a and X? are adjacent cells with common edge 

= X"nW. 

This is the closure of a Jordan open arc C al3 = C a P \ V(X). The mappings 

are ^°°-diffeomorphisms but they do not necessarily match smoothly along 
the edges. We shall now produce a new cell X a/3 , a daughter of the adjacent 
cells X a and X 13 , such that 

C al3 C X afi C x a u c al3 u x p . 
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To construct X a ^ we look at the adjacent dyadic squares T a and T^ 3 in 
Y. The intersection T a n = h{C a P) is a closed interval. Let R be a 
number greater than the length of h(C a ^) to be chosen sufficiently large 
later on. There exist exactly two open disks of radius R for which h(C a ^) is 
a chord. Their intersection, denoted by C a ^, is a symmetric doubly convex 
lens of curvature Thus C a ^ is enclosed between two open circular arcs 

7 a/3 = T a n5£ a/3 c &nd y?a = T /3 n ^£"/3 c T /3. Note that £*P = £0a 

but -y a/3 ^7 /3a . We call 

(3.2) = fr^Z; ' 3 ), a daughter of the adjacent cells 36° and %P . 

As the curvature of the lens £ al3 approaches zero, the area of X afS tends to 
0. This allows us to choose R so that 

( 3 - 3 ) II v Mj^(x-/3) < 5 ; 2a+/3 ■ 

The lenses C a/3 are disjoint because the opening angle of each lens (the 
angle between arcs at their common endpoints) is at most tt/3 and their 
long axes are either parallel or orthogonal, see Figure [3j Therefore, the 
cells X a/3 = f^[ l (C a P) are also disjoint. However, their closures may have a 
common point that lies in V(X). The boundary of X a/3 consists of two open 
arcs 

T a ? = X a ndX a/3 and T? a = X? ndX a ? 

plus their endpoints. These open arcs are ^°°-smooth because they come as 
images of the circular arcs enclosing the lens £ a P under a ^°°-diffeomorphism. 




Remark 3.1. In what follows we shall consider only the pairs (a, /3) of indices 
a = 1, 2, . . . and j3 = 1,2,... which correspond to adjacent cells. Such pairs 
will be designated the symbol aj3. 
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1 o- 



Step 2b. (Replacements in X a/3 ) The lenses C a/3 CY are convex, so with 
the aid of Proposition 2.6 and Theorem 2.5 we may replace hi : X a @ — > C a/3 
with the p-harmonic extension of hi : dX a ^ — > dC a @ . We do this, and denote 
the result by hf : X a)3 -»■ only on the cells in which h x : X^X^UX^ -»■ 
M 2 is not a ^^-diffeomorphism. In other cells we set h 2 = hi- In either 
case hf e h 1 + A {X a ' 3 ) so we define 

h 2 = hi + y^[^ 

Thus we have 

(A 2 ) h 2 - hi e A>(X). 

The advantage of using ^2 in the next step lies in the fact that it is not only 
a ^°°-diffeomorphism in every cell, but also is ^"^-smooth with positive 
Jacobian determinant, up to each edge of the cells created here. These edges 
are ^°°-smooth open arcs. By cells created here we mean not only X a/3 
but also those obtained from the parent cell X a by removing the adjacent 
daughters; that is, 

X a \\JX af} , a = 1,2,... 

aB 

See Figure |4j The estimates of h 2 run as follows. By (3.1) we have, 

e 

5' 



{B 2 



\h 2 



hAW(X) ^ sup{diam£ a/3 } ^ sup{diamY a } ^ 



a 3 



In view of the minimum p-harmonic energy principle, we have 
\\Vh 2 - Vhi\\ C pn 



yjiv^2 - v^iii £ p( U3ea ^) 

aB 



5? 



aB 



V/i2||/;p(x a ' 3 ) "I" 1 1 ^^1 !!/>(:£ 
2f 



2 ^ ^ || £P(3"a P) 

aB 



E 2 " 



a-B 



a/3 



by (3.3). Hence 
(C2) 



\Vh 2 - Vh 



l\\CP(X) 



The minimum energy principle also yields estimate 



In particular 
(D 2 ) 



IVTfc IIP 

I V "-211^(^^/3) 



< llV^lH^p^a/S) + l|V^l||^ 



£p(x\uat a ' 3 ) 



l|Vfti||£ 



Iv^lUprx) ^ I|v^i||£p( 



completing the proof of Step 2. 
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Note that hi is locally bi-Lipschitz in X\ V(X). The exceptional set V(X) 
is discrete. 

Step 3 

We shall now merge all the adjacent cells together, by smoothing hi 
around the edges T a/3 C X a . To achieve proper estimates we need to re- 
move small neighborhoods of all vertices, outside which hi is certainly locally 
bi-Lipschitz. 

Step 3a. First we cover the set C(Y) of corners of dyadic squares by disks 
D c centered at c G C(Y). These disks will be chosen small enough to satisfy 
all the conditions listed below. 

(i) diamO c < e/5 for every c G C(Y), 

(ii) Y, I |Vfor< (|jY, where Ft, = ^' 1 (Dc) J c = h 2 (v) = h(v). 
uev(X) F " ^ ' 

Denote by X Q = X\UF„. We truncate each edge T a ^ near the endpoints 
by setting 

(3.4) = r a/3 nx D . 

These are mutually disjoint open arcs; their closures are isolated continua 
in X\ V(X). This means that there are disjoint neighborhoods of them. We 
are actually interested in neig hborhoods U a/f cF of T° P that are Jordan 
domains in which To C are ^°°-smooth crosscuts with two endpoints 
in &U a P , see Section^ It is geometrically clear that such mutually disjoint 
neighborhoods exist. Now the stage for next substep is established. 

Step 3b. (^^-replacement within U a/3 ) It is at this stage that we will 
improve hi in to a "la 00 -smooth diffeomorphism with no harm to the 



previously established estimates for hi. The tool is Proposition 2.8 As 
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always, we shall make no replacement of 1%2 '■ — > T a if it is already c € co - 
diffeomorphism. Recall that we have a bi-Lipschitz mapping fi2 : U a ^ — > 
fi>2(3i a ) = T a that takes the crosscut C onto a circular arc. Denote 
the components H'f = V af} \ and U ^ = U Q/3 n Furthermore, we 
have 



\D%2\ ^ and det-Dfo ^ 



1 



for some m a p > 



ra/3 



T a and fi 2 : 



T «/3 



T a are 



on each component. The mappings H2 : 
^^-diffeomorphisms up to r" . In accordance with Proposition 2.8 we find 
a constant M Qj g such that: whenever open set V a/3 C U a/3 contains the 
crosscut To^ there exists a homeomorphism : ^(U"^) which 

is a ^°°-diffeomorphism in U°^, with the following properties 

• hf = n 2 on (u^\v^)ur^ ; 

• \Vh% \ sC M aP and det Vt% p > ^ in U a/3 . 

Since M Q( g does not depend on Y a @ it will be advantageous to take neigh- 
borhoods Y a P of r° thin enough to satisfy 



|v a ^| < 1 



m a/ 3+M al 3 



and also \Y a ^\ ^ ^ 



2<*+0M afj ■ 



Note that f£ P ,fh G W l ^{W p ) C ^ 1 > p (U a/3 ) and = h 2 on <9U Q/3 , 



so 



we have 



Step 3c. We now define a homeomorphism ^3 : X Y by the rule 

J a/3 

£2 



m 

inX\ 1JU Q/3 . 

a/3 



Obviously, ^3 is a ^^-diffeomorphism in X and ^3 — 1%2 G /^> 1,P (X ). Since 
^3 coincides with ^2 outside X we have ^3 = ^2 + [^3 — ^2)0- Hence 

(A 3 ) ^3-n 2 G.4o(X). 

Then, for every x G X, 



\h 3 (x) - h 2 (x)\ < 



diamft 2 (U a/3 ) 
0, 



for x G U Q 
otherwise 



< diamT" < 



meaning that 



1^3 - foll^X) 
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The computation of p-norms goes as follows 



\vn 3 - vft 2 ||^ m = E I |v^ 3 - vn 2 \p 

a /3 



^Ei v ° 



) ~t~ || v^H^y*/ 3 ) 



at 



Hence 

(c 3 ) IIv^-vm^*)^- 

In the finite energy case, when || V^|| ^fp(x) < °°> we observe that 

ll V ^3||.£ , P(X\UV«' 3 ) = l|V^2||jfP(X\UV a ' 9 ) ^ ll V M^P(X)- 

Therefore, by triangle inequality, 

l|V^3||jSfP(X) ^ l|V^2||jfP(X) + Ell V ^3||^p(y^/9) 

a/3 

< liv^ll^(x) + El va/3 I • H V M^( V ^) 

a/3 

^l|V^||^(X)+E 2 a +/ 3 M 

a/3 ap 

which yields 

(AO HV^Ibnx) < l|V^||^( X) + <5. 

The third step is completed. 

Step 4 

We have already upgraded the mapping h to a homeomorphism h 3 : X 
Y that is a ^°°-diffeomorphism in X D = X \ Uuev(x) wnere F v are small 
surroundings of the vertices of cells. Their images h 3 (¥ v ) = fi2(¥ v ) = D c are 
small disks centered at c = h(v). In Step 3a, one of the preconditions on 
those disks was that diamD c < e/5. Furthermore, the closed disks D c are 
isolated continua in Y for all c 6 C(Y), so are the sets ¥ v in X. We shall now 
consider slightly larger concentric open disks D' c D D c , c G C(Y), and their 
preimages = 1 (W C ) cX,v = h _1 (c) G V(X). The annulus D' c \ W c will 
be thin enough to ensure that W c are still disjoint, 

diamD' c < ^ for all c G C(Y) 

and 

e \p 



«ev(x) 
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Figure 5. Neighborhoods of vertices. 



Let F' v , v G V(X), denote the boundary of ¥' v . These are ^°°-smooth 
Jordan curves. We now define a homeomorphism : X Y by perform- 
ing p- harmonic replacement of mappings h 3 : ¥' v W c , whenever such a 
mapping fails to be ^°°-diffeomorphism. Thus every h^: ¥' v W c is a 
*^°°-diffeomorphism up to 1^. Moreover fi 4 efi 3 + Wo' p (F' c ), so 



h A -H 3 £i„(X). 



For every x € X, we have 



\h 4 (x) - h 3 (x)\ < 



diam D 



inF;, c = 
otherwise 



h(v) 



Hence 



11^4 - M%?(X) 
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By virtue of the minimum energy principle we compute the p-norms 



— ^3\\%p(x) ~ E 11^4 — ^3 I|P 



Hence 

(C4) 



< E [HM-SfP(F{,) + II^3||^p(F^)] P 

ueV(X) 
ueV(X) 

^ 2 2p [\\^3\\%> P (f> v \ ¥v ^ + \\^3\\%- P (Y v ) 



< 2 2p " 1 

< 2 2p 



«ev(x) 

e \p 
20 



+ E life 



20/ ~~ U> 



1^4 — ^: 



Again by minimum energy principle we find that 

(A) 1 1 ^4 1 1 _^( X ) < IIMjsn>(X)- 

Just as in the previous steps, condition (£4) remains valid, finishing Step 4. 



Step 5 

The final step consists of smoothing in a neighborhood of each smooth 
Jordan curve T' v , v G V(X). We argue in much the same way as in Step 3, 



but this time we appeal to Proposition 2.9 instead of Proposition 2.8 By 
smoothing H4 in a sufficiently thin neighborhood of each T' v we obtain a 
^^-diffeomorphism ^5 : X Y, 



(A 5 



h 5 -h±£ «4o(X). 



(B 5 ) 

(OO 
(A) 



|% - ^4|k(X) ^ 



11^5 - M SfpfX) < 



ll^5||jSfp(X) < ||^4||jSfP(X) + & a 
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4. Open questions 
Question 4.1. Does Theorem 1 1 . 1 1 extend to n = 3? 

Question 4.2. A bi-Sobolev homeomorphism h: X Y is a mapping 
of class W l ' p {X,Y), 1 < p < oo, whose inverse 

h -i . Y ont^ x belongg t() 

a Sobolev class ^ 1,<? (Y,X), 1 ^ g < oo. Can /i be approximated by bi- 
Sobolev diffeomorphisms {/i^} so that hi — > h in ^ 1,P (X, Y) and /i^T 1 — > h^ 1 
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